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A  Finite  Element  Approach  for  Multidimensional  Inverse 

Heat  Conduction 


Edward  B.  Coy1  and  Matthew  Bergkoetter2,  Stephen  A.  Danczyk3,  Edgar  Felix4 
Air  Force  Research  Laboratory,  Edwards  AFB,CA, 93524 


An  efficient  technique  for  mapping  thermal  boundary  conditions  is  described  and 
demonstrated.  The  technique  is  based  on  a  piece-wise  polynomial  approximation  where  the 
Lapiacian  derivatives  in  space  are  constrained  using  the  heat  equation.  Measured  values  for 
the  Lapiacian  are  obtained  from  temperature  rate  measurements  from  sensors  embedded 
within  a  body.  The  technique  has  been  implemented  in  a  digital  signal  processor  and  is  able 
to  provide  real-time  data  on  thermal  boundary  conditions  over  a  surface.  The  technique  is 
adaptable  to  complex  geometry.  In  this  paper  the  technique  is  applied  to  a  study  of  the 
injector-wall  interactions  in  a  laboratory  scale  liquid  rocket  engine. 
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Nomenclature 

=  a  constant  in  a  polynomial,  also  a  dimension 
=  vector  of  data 
=  a  matrix  of  polynomial  terms 
=  a  vector  of  polynomials 
=  heat  flux 
=  temperature 
=  time 

=  plate  thickness 
=  dimensional  coordinates 
=  a  complex  constant 
=  angular  frequency 
=  complex  conjugate 
=  frequency  domain  variable 
=  transpose 
=  inverse 


I.  Introduction 

Local  concentrations  of  heat  flux  and  temperature  are  of  interest  in  aerospace  applications  including  injector-wall 
interactions  in  liquid  rocket  engines,  shock-shock  interactions  in  external  flow  over  hypersonic  vehicles, 
internal  and  external  flow  over  blades  and  vanes  in  gas-turbine  engines,  and  deposition  of  directed  energy.  Accurate 
characterization  of  the  distributions  may  be  required  for  validation  of  a  model,  optimization  of  a  design,  or  control 
of  an  operational  system.  There  are  two  basic  methods  in  use  today.  In  the  calorimetry  approach  a  heat  exchanger  is 
incorporated  into  the  wall  of  the  device  and  the  temperature  rise  of  the  working  fluid  is  used  to  determine  the  local 
heat  flux.  This  method  is  limited  to  applications  where  the  heat  flux  is  slowly  varying  in  time  and  changes  in  one 
direction  only  along  the  surface.  Various  assumptions  and  approximations  are  required  to  infer  the  surface 
temperature;  although,  in  some  cases,  temperature  sensors  can  be  incorporated  to  give  direct  measurements.  To 
obtain  high  spatial  resolution  with  a  calorimeter,  specialized  manufacturing  processes  are  required.  The  other  major 
class  of  methods  is  appropriate  for  transient  situations  which  commonly  occur  in  hypersonic  wind  tunnel  test  articles 
and  laboratory  scale  combustion  devices  where  the  test  articles  are  designed  as  heat-sink  devices.  Surface  junction 
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thermocouples  have  been  developed  which  match  the  thermal  characteristics  of  the  test  articles  and  can  be  treated  as 
non-intrusive  to  the  flow  of  heat.  From  the  temperature  history  a  heat  flux  history  can  be  constructed.  In  some 
situations  the  temperature  sensors  must  be  embedded  in  the  wall  at  some  distance  from  the  surface.  In  this  situation, 
inverse  methods  must  be  used  to  reconstruct  the  surface  heat  flux  and  temperature1.  The  mathematical  methods  for 
performing  the  analysis  for  surface  junction  and  embedded  sensors  have  been  extensively  studied  and  are  highly 
developed;  however  the  most  widely  used  methods  assume  that  the  flow  of  heat  is  one-dimensional.  The  potential 
errors  in  this  approach  can  be  readily  appreciated  when  one  considers  the  situation  that  occurs  when  heat  is 
deposited  by  a  concentrated  source  such  as  a  laser.  When  energy  is  deposited  adjacent  to  the  sensor,  heat  diffuses  in 
the  plane  of  the  surface  causing  a  temperature  rise  at  the  location  of  the  sensor  which  is  recorded  incorrectly  as  heat 
flux  at  that  point.  Due  to  a  plethora  of  such  applications,  research  in  this  area  has  shifted  to  the  development  of 
solution  techniques  for  spatially  multi-dimensional  and  transient  problems. 

The  mathematical  characteristics  of  the  inverse  heat  conduction  problem  are  well  known.  Changes  in 
boundary  conditions  cause  attenuated  and  time-lagged  responses  within  a  body.  In  the  presence  of  noise  and 
measurement  error,  the  problem  becomes  ill-posed  and  solutions  can  be  biased  and  unstable.  Successful  inverse 
techniques  have  been  developed  which  rely  on  the  use  of  information  from  future  time  steps  to  ensure  that  changes 
in  the  boundary  conditions  at  the  current  time  step  have  had  time  to  reach  the  sensor  locations.  Stabilization  has  also 
been  accomplished  through  Tikhonov  regularization.  Other  techniques  utilize  low-order  basis  functions  to  dampen 
instabilities.  Most  techniques  rely  on  an  iterative  solution  wherein  an  initial  estimate  for  the  boundary  condition  is 
made,  followed  by  a  forward  calculation  to  establish  the  effect  at  the  sensor  locations,  and  then  a  refinement  in  the 
boundary  condition  with  repetition  until  a  stopping  criteria  is  satisfied.  The  need  to  solve  the  forward  problem 
multiple  times  at  each  time  step  in  a  transient  analysis  results  in  lengthy  run  times  and  must  be  performed  in  post¬ 
processing.  In  addition  to  the  destabilizing  effects  of  noise  and  measurement  error  and  the  time  required  to  complete 
a  calculation,  practical  methods  must  also  contend  with  the  temperature  dependence  of  thermal  properties,  complex 
geometries,  the  need  to  accurately  model  the  effects  of  temperature  sensor  installations,  have  the  ability  to  handle 
transient  and  steady  state  applications,  and  be  intelligible  to  those  responsible  for  applying  the  technique. 


II.  Theory 


In  the  one-dimensional  approximation2,  the  temperature  profile  within  a  body  can  be  approximated  with  a 
polynomial. 


T(z)~  P„(z)  =  c0  +  c1z  +  ...  +  c„z"  (1) 

In  general,  n+1  measurements  are  required  to  solve  for  the  n+1  coefficients.  Multiple  measurements  can  be 
obtained  from  a  single  sensor  if  we  make  use  of  the  heat  equation  to  relate  spatial  second  derivatives  with  time 
derivatives.  If  we  consider  two  sensors  that  are  placed  at  two  depths,  and  a  P3  approximation,  we  have  the  following 
solutions. 
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Higher  order  approximations  can  be  constructed  using  spatial  fourth  derivatives  and/or  a  larger  number  of 
sensors.  In  the  one-dimensional  approximation  the  spatial  second  derivatives  are  obtained  from  the  heat  equation. 

d2T  _  1  dT 
dz 2  a  dt 

A  correction  for  the  effects  of  temperature  gradients  normal  to  the  z-direction  can  be  made  using  adjacent 
sensors  to  form  a  finite-difference  approximation  to  the  Laplacian  derivative. 
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d2T  _  T{x  -l,z)~  2 T(x,  z)  +  T(x  +  /,  z) 
dx2  ~  l 2 


d2T  '  1  dT  T(x-l,z)-2T(x,z)  +  T(x  +  l,z) 
dz2  a  dt  l 2 


The  accuracy  of  the  correction  for  lateral  diffusion  can  be  evaluated  by  comparing  with  analytical  solutions.  A 
suitable  test  case  is  a  semi-infinite  body  exposed  to  a  step  change  in  heat  flux  over  half  the  surface  and  insulated 
over  the  other  half3.  The  dimensional  quantities  used  in  the  test  case  were  chosen  to  be  representative  of  a  device 
and  test  conditions  described  later  in  this  paper.  The  edge  of  the  heated  region  is  located  at  the  centerline,  x=0.  The 
model  contains  nine  pairs  of  temperature  sensors  at  depths  ofzy=1.65  mm  andz?=5.21  mm,  with  a  lateral  spacing  of 
/=4.23  mm.  The  heat  flux  used  in  the  simulations  was  5  W/mm2.  The  material  properties  of  the  solid  were  reference 
values  for  copper,  £=400  W/m-K  and  a=l  17  mm2/s. 


<1 


Figure  1  Validation  case.  Semi-infinite  body  with  uniform  heat  flux  over  half  the  face. 

Figure  (2)  contains  heat  flux  results  for  the  one-dimensional  approximation,  equation  (4),  and  the  corrected  ID 
approximation,  equation  (6),  which  is  labeled  2D  in  the  figure.  There  is  a  significant  improvement  achieved  by 
accounting  for  the  lateral  diffusion  of  heat.  At  the  centerline  where  the  edge  of  the  heated  region  is  located,  the 
predicted  heat  flux  is  the  average  value.  Figure  (3)  shows  the  prediction  of  surface  temperature.  The  2D 
approximation  achieves  a  better  agreement  with  the  analytical  solution,  although  surface  temperature  is  a  less 
sensitive  quantity  to  predict. 
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Figure  2  Heat  Flux  t=l  s 
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Figure  3  Surface  Temperature  t=l  s 


Figures  (4)  and  (5)  show  simulated  effects  of  noise  in  the  temperature  measurements  and  the  sensor  locations. 
The  random  error  in  temperature  was  assumed  to  be  1  K  and  the  random  error  in  sensor  locations  in  both  x  and  z 
directions  was  assumed  to  be  0.1  mm.  These  are  representative  values  based  on  experience.  Twenty  simulations 
were  conducted  with  random  errors  generated  from  a  uniform  distribution.  The  estimated  error  in  the  heat  flux  for 
the  simulated  conditions  is  in  the  range  of  +/-  0.5  W/mm2  and  the  estimated  error  in  the  surface  temperature  is  +/-  5 
K.  These  values  depend  on  many  factors.  To  get  as  accurate  an  estimate  as  possible  the  simulated  conditions  need  to 
mimic  the  test  conditions.  In  general,  the  fractional  uncertainty  decreases  with  increasing  heat  flux. 


Figure  4  Noised  data  <7v=0.1mm,  aT= IK 
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Figure  5  Noised  data,  <rv=0.1mm,  <7r=lK 


III.  Validation 

Substitution  calibration  of  heat  flux  sensors  for  liquid  rocket  engine  applications  is  not  possible  because  there  is 
no  reference  standard  available.  Estimates  of  accuracy  must  rely  on  analysis  of  the  sources  of  error.  In  addition  to 
the  factors  already  discussed  in  Section  II,  the  sources  of  error  include  systematic  errors  in  the  temperature 
measurements,  the  thermal  transport  property  correlations,  and  the  assumptions  used  in  deriving  the  mathematical 
model.  In  this  section  we  describe  tests  we  used  to  quantify  these  sources  of  error.  The  method  is  based  on  the 
following  analytical  solution  to  the  one-dimensional  heat  equation  for  a  semi-infinite  domain  with  a  steady-periodic 
heat  flux  at  the  surface. 


T(x,t)  =  T0+^- 
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If  we  can  show  our  assumptions  are  valid  for  the  conditions  of  equation  (7),  then  they  must  also  be  valid  for 
equations  (2)  and  (3).  The  advantage  of  using  equation  (7)  is  that  we  can  avoid  the  measurement  of  an  absolute  level 
of  heat  flux  and  instead  base  the  validation  on  a  measurement  of  thermal  diffusivity  which  is  well  known.  The 
measurement  is  made  by  exposing  a  heat  flux  sensor  to  a  steady  periodic  heat  source  and  feeding  the  output  of  the 
temperature  sensors  into  a  lock-in  amplifier  capable  of  measuring  phase  angle.  The  thermal  diffusivity  is  obtained 
from  the  following  expression  derived  from  the  phase  angle. 
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If  the  measured  thermal  diffusivity  does  not  agree  with  the  reference  value,  the  errors  may  be  due  to 
thermocouple  placement,  thermocouple  time  response,  multidimensional  heat  flow,  or  impurities  in  the  material. 
Equation  (8)  will  not  identify  the  source  of  error  but  it  is  unlikely  that  there  will  be  cancelling  errors  that  lead 
fortuitously  to  a  correct  thermal  value  of  thermal  diffusivity. 

In  the  following  figure  we  show  an  example  of  errors  due  to  multi-dimensional  heat  flow.  A  carbon-tube  furnace 
(Thermogauge  Instruments  Model  HT-7000A)  was  used  as  a  heat  source  and  an  optical  chopper  (New  Focus  3501) 
was  used  to  create  a  square  wave  input  at  0.5  Hz.  The  phase  angles  were  measured  with  a  digital  lock-in  amplifier 
(Signal  Recovery  7265).  The  sensor  element  was  a  cylindrical  rod  of  copper  (Alloy  101),  12.5  mm  in  diameter  and 
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50  mm  long,  pressed  into  a  copper  cylinder,  19.1  mm  in  diameter  and  50  mm  long,  with  two  type-K  thermocouples 
pressed  into  grooves  on  opposing  sides  of  the  inner  rod.  The  sheaths  were  stripped  back  exposing  the  thermocouple 
wires.  An  alumina  insulator  was  placed  over  the  wires  to  avoid  contact  between  the  wires  and  the  sheath,  between 
the  wires  themselves,  or  between  the  wires  and  the  base  material  at  an  incorrect  point.  The  wires  were  pressed  into 
circumferential  grooves  in  the  inner  rod  and  the  junctions  were  formed  through  the  base  material  in  the  intrinsic 
configuration.  As  the  opening  in  the  chopper  wheel  swept  across  the  surface  of  the  sensor,  it  first  passed  over  one 
thermocouple  location  and  then  the  other,  thereby  creating  a  flow  of  heat  that  was  not  one-dimensional.  Rotating  the 
sensor  on  its  axis  caused  the  relative  position  of  the  two  sensors  to  shift  and  resulted  in  a  varying  value  for  the 
measured  thermal  diffusivity. 


Figure  4  Erroneous  thermal  diffusivity  due  to  multidimensional  heat  flow 


The  experiment  was  repeated  with  the  optical  chopper  replaced  with  a  fast  shutter  (Uniblitz  VMM-T1  driver  and 
CS45  shutter).  In  this  case,  both  thermocouple  locations  received  simultaneous  exposure  and  the  thermal  diffusivity 
results  were  much  improved.  These  results  validate  the  method  with  respect  to  the  method  used  to  install 
thermocouples,  the  thermal  property  correlations,  and  the  assumptions  used  in  the  derivation  of  the  one-dimensional 
model.  The  results  also  illustrate  the  importance  of  accounting  for  multidimensional  heat  flow.  The  results  of 
Section  II  above  validate  the  extension  to  multidimensional  heat  flow. 
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Figure  5  Results  with  one-dimensional  heat  flow. 
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III.  Implementation 


The  code  that  implements  equations  (2)  and  (3)  in  the  test  facility  runs  on  a  Pacific  Instruments  model  6042  card, 
which  is  built  around  a  Texas  Instruments  Model  TMS320C6711  Floating  Point  Digital  Signal  Processor  capable  of 
900  million  floating  point  operations  per  second  at  a  clock  rate  of  150MHz.  The  PI  card’s  firmware  includes  a  finite 
impulse  response  filter  function  as  well  as  basic  matrix/array  addition  and  multiplication  which  are  accessible 
through  the  PI660  GUI  software  from  Pacific  Instruments.  For  the  sake  of  convenience,  the  current  implementation 
is  built  from  these  basic  operations  in  the  PI660  environment,  although  a  lower-level  implementation  compiled 
using  the  TI  processor’s  software  development  kit  (SDK)  could  be  expected  to  improve  efficiency.  The  current 
method  has  been  demonstrated  with  32  channels  at  a  rate  of  1000  Hz.  The  major  steps  in  the  implementation  are  as 
follows: 

1 .  Load  pre-calculated  constants  before  the  system  starts  taking  data 

2.  Load  data  from  DAQ  backplane  to  DSP  memory 

3.  Calculate  temperature  values  from  raw  voltage 

4.  FIR  filter  data  (smoothing  and  differentiating) 

5.  Estimate  thermal  diffusivity  a  at  each  sensor  location  based  on  linear  fitting  coefficients  for  a(T) 

6.  Matrix  multiplication  for  T(0)  and  quantity  in  parenthesis  in  eq.  (3) 

7.  Estimate  thermal  conductivity  k  for  each  of  the  surface  temperature  values  using  a  linear  fit  of  k(T) 

8.  Find  q  on  the  surface  using  k(T(0)) 

9.  Send  results  back  to  DAQ  system 


IV.  Application 

The  target  application  for  this  effort  was  the  characterization  of  injector-wall  interactions  in  liquid  rocket  engine 
combustion  chambers.  Most  chamber  failures  occur  due  to  localized  concentrations  of  heat  flux.  High  heat  flux  and 
high  oxygen  concentrations  impinging  on  the  chamber  wall  can  lead  to  a  phenomena  known  as  blanching4.  We  are 
interested  in  creating  a  capability  for  screening  candidate  injector  designs  for  this  type  of  behavior  in  lab-scale  tests 
to  avoid  costly  redesigns  of  full-scale  engines. 

For  this  study  we  screened  two  injector  types.  The  first  was  a  gas-centered  swirl  co-axial  element  (GCSC)5'6 
operating  with  gaseous  oxygen  and  RP-2,  a  liquid  hydrocarbon  propellant.  The  second  was  an  array  of  25  shear,  co¬ 
axial  elements  operating  with  gaseous  oxygen  and  gaseous  hydrogen.  In  both  cases  we  characterized  the  heat  flux 
patterns  produced  by  the  injectors  in  an  as-built  condition  and  after  “defects”  were  deliberately  created.  The  gas- 
centered  swirl  coaxial  element  design  is  shown  in  Figure  (7).  The  cylindrical  passage  is  referred  to  as  the  post  and  is 
for  the  gaseous  oxygen.  The  post  includes  three  liquid  injection  orifices  near  the  downstream  end.  The  liquid  enters 
the  central  cylinder  tangentially  to  generate  a  swirling  liquid  film  around  the  inner  diameter  of  the  post.  The  liquid 
film  is  thus  subjected  to  a  combination  of  cross-flow  shear  and  centrifugal  forces.  The  liquid  is  stripped  from  the 
film  inside  the  element  by  the  central  gas  jet,  which  entrains  the  droplets  and  transports  the  resultant  spray 
downstream.  In  this  study  we  blocked  one  of  the  three  fuel  injection  orifices  in  an  attempt  to  create  a  non-uniform 
spray  pattern  and  wall  heat  flux. 
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Figure  7  Lab-scale  gaseous-oxygen,  co-axial  swirled  injector 

The  25  element  array  of  shear  co-axial  injector  elements  is  shown  in  Figure  (8).  The  gaseous  oxygen  enters  from 
the  left.  It  first  encounters  a  distribution  plate  and  then  enters  the  oxygen  posts.  Each  oxygen  post  contains  a 
precision  orifice  pressed  into  the  inlet.  The  hydrogen  enters  through  the  openings  in  the  top  and  bottom  and  flows 
through  passages  in  the  oxygen-post  support  plate  and  then  the  annular  gaps  between  the  oxygen  posts  and  the  face 
plate.  To  create  non-uniformity,  the  inlet  orifice  in  the  oxygen  post  closest  to  the  wall  between  the  centerline  and 
wall  was  blocked. 


Figure  6  Lab-scale  shear  co-axial  injector  elements 

The  measurement  device  was  located  38  mm  from  the  face  plate.  It  consisted  of  a  block  of  copper  with  nine  pairs 
of  embedded  thermocouples.  The  pairs  were  oriented  perpendicular  to  the  direction  of  gas  flow  to  capture  streaks. 
The  thermocouple  junctions  were  formed  through  the  copper  base  material  using  the  technique  of  intrinsic 
thermocouples  to  ensure  a  fast  time  response.  An  alumina  insulator  was  used  to  ensure  that  the  only  point  of  contact 
for  the  thermocouple  wires  was  on  the  copper  surface  at  the  bottom  of  the  hole.  The  thermocouple  wires  were 
terminated  on  an  Iso-Tech  TRU  100  Multi-Channel  TC  Referencing  Unit  and  data  were  recorded  on  Pacific 
Instruments  Model  6013  cards  in  a  PI  660  DAQ  system  at  1000  Hz  with  16  bit  resolution  and  conditioned  with  25 
Hz  4-pole  Bessel  filters. 

The  heat  flux  and  temperature  distributions  for  the  two  injectors  in  the  as-built  and  defect  configurations  are 
shown  in  Figures  (7)-(  1 0).  The  data  shown  was  taken  immediately  prior  to  the  end  of  the  run.  The  heat  flux  profile 
for  GCSC  is  highly  uniform  for  the  as-built  configuration  and  skewed  for  the  defect  configuration.  The  co-ax  is 
skewed  for  both  configurations.  Surface  temperatures  follow  the  same  trends  but  are  skewed  for  all  cases  due  to 
residual  effects  of  the  torch  igniter  that  created  an  initially  skewed  temperature  profiles  as  it  fired  across  the  face  of 
the  injector  during  startup.  There  are  no  significant  streaking  patterns  observed  under  any  conditions.  In  future  work 
we  will  compare  these  measurements  with  predictions  using  a  CFD  code. 
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Figure  7  Heat  flux  distributions  for  GCSC  injector 


-♦-Co-ax  as-built  «**Co-ax+defect 


Figure  8  Heat  flux  distributions  for  shear  co-ax  injectors 
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-♦-GCSC  as-built  -»-GCSC  +  defect 


Sensor# 

Figure  9  Surface  temperature  distributions  for  GCSC  injector 


-♦-Co-ax  as-built  -«-Co-ax+defect 


Figure  10  Surface  temperature  distributions  for  shear  co-ax  injectors 


V.  Multidimensional  Approach 

In  this  final  section  we  provide  some  preliminary  results  on  an  extension  of  the  general  method  described  in 
Section  II  to  higher  dimensions.  We  utilize  the  strategy  of  domain  decomposition.  Finite  element  methods  are  highly 
developed  for  forward  problems.  They  are  used  because  the  basic  approach  of  subdividing  the  domain  into  cells  and 
approximating  the  solution  over  each  cell  with  a  polynomial  is  readily  adaptable  to  complex  geometries  and  is 
mathematically  tractable.  However,  there  are  some  distinct  differences  between  forward  and  inverse  problems  that 
result  in  some  changes  in  the  basic  strategy.  In  solving  forward  problems  it  is  possible  to  generate  meshes  that  are 
limited  in  size  only  by  the  memory  capacity  of  the  computer.  It  is  possible  to  finely  discretize  regions  where  the 
function  behavior  is  complex  such  that,  within  each  individual  cell,  the  function  can  be  approximated  adequately 
with  a  simple  linear  basis  function.  In  inverse  problems,  as  developed  here,  the  grid  corresponds  to  a  network  of 
sensors.  In  practical  applications,  grids  with  tens  or  hundreds  of  nodes  are  possible,  but  probably  not  thousands,  and 
certainly  not  the  millions  that  are  now  common  in  forward  problems.  This  provides  a  strong  impetus  for  developing 
schemes  that  use  high  order  basis  functions.  Coy  showed  how  basis  functions  of  arbitrary  order  can  be  developed  in 
one  spatial  dimension  using  two  nodes  by  constraining  the  approximating  polynomials  using  the  even  numbered 
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spatial  derivatives  from  the  heat  equation.  Adding  second  derivatives  improved  the  frequency  response  by  three 
orders  of  magnitude.  The  model  could  be  formulated  as  a  digital  filter  and  calculated  in  real  time.  In  this  paper  we 
extend  the  approach  to  three  dimensions. 

To  fix  ideas  we  consider  a  set  of  sensors  embedded  in  a  solid  material  with  smoothly  varying  thermal 
transport  properties.  We  are  primarily  interested  in  constructing  approximations  for  the  temperature  and  heat  flux 
distribution  at  the  surface.  The  temperature  distribution  at  any  point  in  time  within  each  finite  element  is 
approximated  by  a  polynomial. 


Or  simply  as, 


T(x,y,  z,t)  aD  +  QjX  +  a^y  +  a^z  +a4xy  +  asxz  +  a6yz  +  — 


Where, 


rOc,t)  ftf  PT  a 

PT  =  [1  x  y  z  xy  xz  yz  ...] 


(1) 

(2) 

(3) 


a  —  [o0  aL  Qj  Og  Qj  Qg  ...] 


(4) 


A  grid  of  sensors  may  enable  a  finite-difference  estimate  of  the  Laplacian  in  the  plane  of  the  surface.  In  this  case, 
the  second  derivative  normal  to  the  surface  can  be  estimated  as, 


d2T 

dz1 


1  dT 
a  dt 


V\,T 


The  time  derivatives  can  be  obtained  from  the  temperature  records  at  the  measurement  points.  The  set  of  equations 
to  be  solved  can  be  written. 


Or  simply  as, 
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(7) 

(8) 


If  the  number  of  equations  exceeds  the  number  of  terms  in  the  polynomial,  we  can  solve  as  a  least-squares  problem. 

a  =  (GTGy1GTd  (9 

Note  that  i'  GrC)-1Gr  depends  only  on  the  location  of  the  sensors. 

The  performance  of  the  method  can  be  evaluated  by  comparing  with  analytical  solutions.  Cole  has 
developed  Green’s  function  solutions  for  steady-periodic  heating  that  are  suitable  for  this  purpose.  Consider  a 
uniform  slab  of  depth  W  and  infinite  width  with  a  uniform,  single-frequency,  steady-periodic  heat  flux  on  the 
surface  at  y=W,  -a<x<a. 
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The  temperature  inside  the  slab  can  be  expressed  as, 

T  (x,  y,  t)  =  Rea  l  [fix,  y.  &>)  eJ“r  ] 


(10) 
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We  will  use  the  analytical  function  to  obtain  inputs  to  the  finite  element  model,  and  again  to  compare  the 
amplitude  and  phase  of  the  analytical  function  with  predictions  of  the  finite  element  model  at  points  distant  from  the 
sensor  locations,  most  importantly  at  the  surface.  The  model  coefficients  are  real  constants  that  multiply  the  data 
vector  so  we  can  write, 


a  =  ( gtg)~1gt3 

'  5  ' 

t 

a  = 

jvfja  j. 

■}<*>Tjaw 

The  complex  temperature  predicted  by  the  finite  element  model  can  be  written. 


(11) 
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The  complex  heat  flux  in  the  z-direction  is, 


#  v  dP‘ 
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02 

The  gain  and  phase-lag  of  the  finite  element  model  relative  to  the  analytical  solutions  are  the  following, 
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At  locations  .r  =  +a  there  are  gradients  of  temperature  in  the  x-direction  that  we  would  like  to  capture.  We 
take  as  a  length  scale  the  distance  from  the  heated  surface  to  the  first  sensor,  Sj.  Thus  we  have, 
u~  =  1.  iSr  =  u;  u. ,  c*  =  a  J-  ,  =  l'/uj,  =  ;r,  J: . x+  =  x  J.  ,z+=z  J:.  Temperature  and  frequency 

are  non-dimensionalized  as  follows:  T+  =  Tk/qc  a,  and  a»+  =  aiSl/a.  We  consider  a  case  where  a+  and  are 
large  and  look  at  the  effect  of  frequency,  a+,  on  the  surface  temperature,  T+ ,  amplitude  and  phase  for  several  values 
of  the  sensor  spacing,  i+. 

Figure  11  illustrates  the  methods  ability  to  reconstruct  a  complex  surface  temperature  distribution.  The  figure 
shown  is  taken  from  a  transient  simulation  that  illustrates  the  methods  ability  to  reproduce  full  images  of  the  thermal 
boundary  conditions  in  real  time. 


Figure  11  Full  panel  simulation  with  three  band  heat  sources  and  piecewise  polynomial  reconstruction.  The 
left  panel  is  an  8x8  cell  reconstruction  and  the  right  cell  is  the  analytical  solution. 
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VI.  Conclusion 


Direct  measurement  of  complex,  transient,  thermal  boundary  conditions  is  often  not  possible  in  aerospace 
applications  and  inverse  techniques  using  embedded  sensors  must  be  used.  Existing  techniques  for  resolving  the 
inverse  problem  are  computationally  intensive  and  cannot  be  performed  in  real  time  or  even  on  a  post-processing 
time  scale  that  supports  test  schedules.  In  this  paper  we  describe  a  technique  that  addresses  all  of  these  technical 
challenges  and  illustrate  its  application  in  the  characterization  of  injector- wall  interactions  in  a  lab-scale  liquid 
rocket  engine  combustion  chamber.  Data  was  acquired  from  two  injector  types  in  an  as-built  configuration  and  with 
defects  introduced  in  an  initial  attempt  to  assess  the  sensitivity  of  the  designs  to  manufacturing  flaws.  In  addition  to 
the  analytical  model,  we  include  an  experimental  characterization  of  the  method’s  accuracy  and  a  description  of  its 
implementation  using  a  digital  signal  processor.  A  further  generalization  and  extension  of  the  technique  to  multi¬ 
dimensional,  transient  problems  for  arbitrary  geometries  is  described. 
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